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Abstract 

This paper gives a definition of g-harmonic functions and shows the 
relation between the g-harmonic functions and g- martingales. It's direct 
to construct such relation under smooth case, but for continuous case 
we need the theory of viscosity solution. The results show that under 
the nonlinear expectation mechanism, we also can get the similar relation 
between harmonic functions and martingales. Finally, we will give a result 
about the strict converse problem of mean value property of g-harmonic 
functions. 
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1 Introduction and Preliminary 

Harmonic function (Am = 0) has a probabilistic interpretation as that if Au = 
on R n , then it(-Bf) is a martingale for any x £ R n (see for example [5]). This 
relation between martingale and harmonic function connects probability with 
potential analysis. It helps us to give probabilistic characterizations for har- 
monic function and more generalized X-harmonic function(see 6,). In 1997, 
Peng([9 ) introduced the notions of g-expectation and conditional g-expectation 
via a backward stochastic differential equation(BSDE) with a generator function 
g. Further, Peng([TO]) introduced the notion of g-martingale which provided an 
heuristic tool to characterize some kind of harmonic function described by ellip- 
tic operator with a nonlinear term g by which we will define g-harmonic function. 
In this paper, with the help of nonlinear Feynman-Kac formula established from 
BSDE (see for example [H]), we will give a probabilistic characterization of the 
g-harmonic function. 

Now we state our problem in detail. Let (f2, J 7 , P) be a probability space en- 
dowed with the natural filtration {Tt}t>o generated by an n-dimensional Brow- 
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nian motion {Bf }t>o- i- c - 

F t = a{B s :s<t}. 

Then we can define a g-martingale by an J-t-adapted process {yt}t>o which 
satisfies the following BSDE for any < s < t: 

y s = yt + J g(y r ,z r )dr - J z r dB r . (1) 
Here g : R x R n — > R, satisfies the condition: 

(HI). g(y, 0) = and the Lipschitz condition: 3C > 0, for any (yi, zi), (y 2 , z 2 ) S 
R x R n we have 

\g(yi, z i) - g(v2,z 2 )\ < C(|yi - 2/2 1 + ki - ^D- 

And the equality ((T|) also can be formulated simply as|llj: 

£ 9 s,t{yt) ■= y s - 

Then we can also get the definition of g-super (sub) martingale when 

£Uvt) < (>) vs. 

This definition derives from the definition of g-expectation in the beginning 
paper Peng[9]. When g(y,z) = the g-expectation is actually the classical 
expectation. Except that g-expectation is nonlinear in general, it holds many 
other important properties as its classical counterpart [5] [J] [TU] [H] . 
Given an n-dimensional Ito diffusion process{Xf } t >o- 

dXf = b(Xf)dt + a(Xf)dB t , (2) 
Xj$ = xeR n , 

where b(x) : R n — > R n , a(x) : R n — > ]J" X " satisfy the Lipschitz condition: 
3C > s.t. 

\b(xi) - b{x 2 )\ + \<y{xi) - <r{x 2 )\ < C\xi - x 2 \ Vxi,x 2 e R n , 

our problem is that: what kind of function u(x) : R n — > R satisfies that u{Xf) 
is a g-martingale for any x G R n l 

This problem also has its classical counterpart: 

First if {Xf} is just the Brownian motion {Bf}, then we have the result 
that when u(x) is harmonic on R n i.e. 

d 2 u 
dxf 

the process u(B£) is a martingale for any x. And conversely if u(x) satisfies that 
u(Bf) is a martingale for any x, then u(x) must be harmonic on R n . The proof 
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may have many editions, here we can give a sketch of one which may induce the 
extension to g-martingale case. 

If u(x) is harmonic on R n , then we use Ito formula to u(B^) and get 

f^l I 1 Pi 1 7 / 

■i i % 

i 

Then we get u(Bf) is a martingale for any x € i2". Conversely if is 
continuous on R n and for any x € i?™, u{Bf) is a martingale, then we have 
£J[it(.B;J:)] = u(a;) for any stopping time r. Particularly for any sphere S{x, r) — 
{y G R n : \y — x\ < r}, we have 

u(x) = E[u(B* )} = [ u(y)da v 

JdS(x,r) 

where T S ^ x r ^ is the exit time of {Bf} from the sphere S(x,r), i.e. 

Ts(x,r) = inf{t > : \Bf - x\ > r}, 

and Uy is the harmonic measure on the dS(x, r). Then from the familiar converse 
of the mean value property for harmonic function, we can get u{x) must be 
harmonic function. 

Further we can extend the Brownian motion {-Bf } to the general diffusion 
process {Xf }: 

If u(x) 6 Cg(i? n ) and satisfies 

then we have u(Xf) is a martingale for any x. The proof also uses the Ito 
formula. But conversely if u{Xf ) is a martingale for any x, we can't conclude 
that u(x) is smooth. Then with additional assumption u(x) G Co(i2") we can 
get that u(x) satisfies the PDE©(see [§]). 

Then naturally we will ask that what happens when we substitute the ex- 
pectation mechanism by the g-expectation mechanism. First we will define the 
infinitesimal generator: 

Definition 1. Let 

£f. t \f(Xf)] -fix) 
Aff(x):=ma y J n \ (4) 

then we call the infinitesimal generator of a diffusion process {Xf } under 
g- expectations. 
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Thanks to the celebrating nonlinear Feynman-Kac formula [5], we can get 
the explicit form of Ag when / 6 C$(R n ) by considering the following type of 
quasilinear parabolic PDE: 

§S - Cu(t, x) - g(u(t, x),u x (t, x)a(x)) = 0, 
u(0, as) = /(a;). 

Where 

Cu(t,x) = J2 bi d^ ^ x ^ + 2^ aaT ^dx~dx~ (6) 

When f eC$(R n ), we assert that 

u(t,x) = £l t [f(X?)] (7) 

is the solution of PDE©. Then under the case t = 0, we get 

Aff(x) = Cf(x)+g(f(x)J x (x)a(x)). (8) 

Then we finish the preliminary and we can introduce our main results. In 
section 2, we give a characterization of g-harmonic function under smooth case. 
In section 3, we characterize it under continuous case, where the differential 
operator is interpreted as viscosity solution. In section 4, we will investigate the 
strict converse problem of mean value property of g-harmonic function evoked 
by its classical counterpart [7]. 



2 Smooth Case 

The equality (jHJ has the implication about the relation between the g-martingales 
and the g-harmonic functions when / e Cq(R 71 ). In fact, the left side of (JSj) is 
related to a g-martingale and the right side is related to a harmonic PDE. At 
first we will give the definition of g-harmonic functions: 

Definition 2. Let f G Cq (i?™). We call it a g- (super)harmonic function w.r.t. 
{X£ } if it satisfies 

Agf(x)(<) = for any x £ R n . (9) 

Then we suffice to construct the relation between the g-supermartingales 
and the g-superharmonic functions. 

Theorem 1. If f(x) G Cq(R u ), then the following assertions are equivalent: 

(1) f(x) is a g-superharmonic function. 

(2) {f(X^)} is a g-supermartingale for any x 6 R n . 
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Proof, (i) (1) =► (2): 

/ G C 2 (R n ), by Ito formula, we can get f{Xf) is still an ltd diffusion process: 
f(Xf)= f{Xf)+ f £f(X?)dr+ I f x (X?)o-{X?)dB r 0<s<t. 

J s J s 

and then we insert the term g(f(X^), f x a(X£)) and get: 

f{Xf) = /(Xf ) - f Cf(X?)dr - f f x a(K)dB r 



f(xn+J g(f(X?)J x a(X r x ))dr- J f x a{X*)dB r 
\c/(X?) + g(f(X?), fMK))}dr 
fix) is a g-superharmonic function, so: 

Cf(X?)+g(f(X?),f x o-(X?)) = Ajf(X*) < 0. 

And then according to the comparison theory of BSDE(cf.[in]), we can get 
{f{Xf)} is a g-supermartingale. 

(ii) (2) => (1): 
By the definition of the Af : 

e§t\f(xf)]-f(x) 

A *w=^ — V — • 

{f{X'f)} is a g-supermartingale, so: 

£° >t [f(X?)]-f(x)<0, 

then 

Aff(x) < 0. 

So we get fix) is a g-superharmonic function. □ 



3 Continuous Case 

If we generalize the requirement of function fix) to be only continuous on R n , 
how we get a function / which satisfies that f{Xf) is a g- martingale for any 
x 6 R n ? With the help of viscosity solution(cf. [3]) we can also refer to the 
quasi- linear second order PDEs. Here we need a lemma due to Peng[8]. 

Lemma 1. Let < t < T and 

u{t,x) = ei T _ t [f{x%_ t )\. 
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Then u(t,x) is the viscosity solution of the following PDE on (0,T) x R n : 

(10) 



§f + Cu(t, x) + g(u(t, x),u x (t, x)a(xj) = 
u(T,x) = f(x) 



Hereg(y,z) and f(x) satisfy: 

(H2). Let F(u,p) = g(u,pa(x)), then 3C > s.t. 

\F(u,p)\<C(l + \u\ + \p\); 
\D u F(u,p)\,\D p F(u,p)\<C, 

and (H3). f(x) is a continuous function with a polynomial growth at infinity. 

Definition 3. Let u(t,x) G C(R x R n ). u(t,x) is said to be a viscosity super- 
solution (resp. sub-solution) of the following PDEUD): 

— + Cu(t, x) + g(u(t, x),u x (t, x)a(x)) = 0, (11) 

if for any (t,x) G R x R n and ip G C 1,2 (R x R n ) such that f(t,x) — u(t,x) and 
{t,x) is a maximum (resp. minimum) point of if — u, 

— (t,x) +&p(t,x) + g(ip(t,x),(p x (t,x)o-(x)) < 0. 
ot 

(resp. ^-(t,x) + Cip(t,x) + g((p(t,x),f x (t,x)o-(x)) > 0.) 

u(t,x) is said to be a viscosity solution of PDE[Tl}) if it is both a viscosity super- 
and sub-solution of ill} ). 

We also consider the viscosity solution of the following type of quasilinear 
elliptic PDE(PT): 

Cu(x) + g(u(x),u x {x)a(x)) = 0. (12) 

We can directly get an relation between the two solutions of (fTTj) and (TT2")) : 

Lemma 2. Let u(t, x) — u(x) for all (t, x) G R x R n , then we have: 

u(t,x) is the viscosity super- (sub-) solution of PDE \11\) <^> u(x) is the viscosity 

super- (sub-) solution of PDE \1°2\) . 

Proof. We only prove the case of viscosity super-solution. The "sub-" case is 
an immediate conclusion of the "super-" case, 
(i). »=}.»: 

For any (io,a;o) G R x R n , and a function (p(x) G C 2 (R n ) which satisfies 
f(x) < u(x),(p(xo) — u(xo), we define (f>(t,x) = ip(x) for all (t,x) G R x R n . 
Then 

— = 0, if(t ,x ) = u(t ,x ), <p(t,x) < u(t,x), 
ot 
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and due to the assumption that u(t, x) is the viscosity super-solution of PDE ([TTj) . 
we have 

— (*o,xo) + Ctf(t ,x ) + g(ip(to,xo),<Px(to,xo)a(x )) < 0, 
at 

i.e. 

£tp(x ) + g(ip(x ), (p x (x Q )a(xo)) < 0. 

So u(x) is the viscosity super-solution of PDE ([T2]l . 
(ii). » 

For any (to,Xo) G i? x H n , and a function tp(t,x) G C 2 (i? x R n ) which 
satisfies (p(t, x) < a;), <^(to; #o) = w(to, a^o)j then 

^(«o,a;o) = 0, (13) 



and due to the assumption that u(x) is the viscosity super-solution of PDE(fT2|). 
we have 

C(p(t ,x Q ) + g{(f{t , x ),<p x (to,x a )a(x )) < 0. 
Combined with (|13p. we get 

— (t ,x ) + £<p(to,x ) + g(tp(t ,xo),(p x (t ,xo)cr(xo)) < 0. 

So u(t, x) is the viscosity super-solution of PDE (fTTj) . □ 

Then we can introduce our main result of this section: 
Theorem 2. We have the following two consequences: 

(i) . For any f(x) G C(R n ), and g(y,z) satisfying (HI), if Van £ i? n , f{Xf) is a 
g-supermartingale, then f(x) is a viscosity super- solution of PDE U2\) . 

(ii) . For any f(x) satisfying (H3), and g(y, z) satisfying (HI) (H2), let f(x) 
is a viscosity super-solution of PDE {Up, then {f(X^)} is a g-supermartingale 
for all x G R n . 

Actually, the consequence (ii) is the answer of our main problem and the 
consequence (i) is the converse of it. But (i) is easier to be proved, so we are 
going to prove (i) at first: 

Proof, (i). 

For any x G R n , let <p G C 2 (R n ), <f(x) = f(x) where x is a maximum point 
of tp — f. It means Vx G ii™, we have <p(x) < f(x). Then from ([5]), we get 

Cip(x) + g((f(x), ip x (x)a(x)) = Afcp(x) 

= 1 . m £?Mx?)\- l p(x) 

t|o t 

£?[<p{X?)]-f{x) 



= lim 



tj.o t 
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According to the comparison theory of BSDE, we get 

£2[<p(Xf)]<£°\f(X*)], 
and with the assumption {/(Xf )} is a g-supermartingale, we can get: 
£ B t [<p{Xf)\ - f(x) < £?[f(X?)} - f(x) < 0. 

Then 



g y(g)=]to gfrm]-/(») 

9 rv ' no t 



i.e. 

£<^(:e) + g((p(x),(p x (x)a(x)) < 0. 
By definition, it means f(x) is a viscosity super-solution of PDE (TT2]) , 
(")• 

We want to prove {/(Xf )} is a g-supermartingale for any a; £ R n . It means 
that we need to prove \/x £ R n and V0 < s < t, we have 

£?, t [/(*?)]< 

Under the assumption, in fact b(x), cr(:r) and z) are all independent of time 
t, so we can get the Markovian property of £ 9 S t , i.e. 

£l t \f(Xt)]=£L.\fW-.)]\v=x.. 

Then we get an equivalence relation: 

{f(Xf)}is a g- (super) martingale for any x e R n 

£°[f(X?)} = (<)f(x) for any t > and x <E R n . (14) 

So we suffice to prove the latter assertion. 

For any T > 0, the assumption f(x) is a viscosity super-solution of PDE tfT^t 
implies that f(t, x) := f(x) is a viscosity super-solution to the following PDE: 

I %l+£u(t,x)+g(u(t,x),u x (t,x)a(x))=Q, 
\ u(T,x) = f(x), 

according to lemma [5] And with the help of lemma [TJ 

u(t,x)=£ 9 0tT _ t [f(X£_ t )} 

is actually the viscosity solution of PDE (TT5)) . Moreover by the maximum prin- 
ciple of the viscosity solution ( see pQ), we can get: 

u(t, x) < f(t, x) for any < t < T. 

Especially, we have 

u(0,x)<f(0,x), 

£ 3 T [f{x¥)] < fix). 

□ 
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Corollary 1. (i). For any f(x) G C{R n ), and g(y,z) satisfying (HI), if\/x g 
R n , f(Xt) is a g-martingale, then f(x) is a viscosity solution of PDEFHfy. 
(ii). For any f(x) satisfying (H3), and g(y,z) satisfying (HI) (HE), let f{x) is 
a viscosity solution of PDE\12\). then {f{Xf )} is a g-martingale for all x £ R n . 

It is an immediate consequence from the theorem [5J 

4 Strict Converse of Mean Value Property 

For classical harmonic function, many generalized results of the converse prob- 
lem of mean value property have been investigated (cf. [5][7])- In [7], 0ksendal 
and Stroock give a technique to solve a strict converse of the mean value prop- 
erty for harmonic functions. Now we will generalize it to the case of g-harmonic 
function. Here the strictness means that for each x £ R n we don't need justify 
that for any stopping time r whether £q t (/(X^)) equals f(x). We only need 
to justify one appropriate stopping time of each x. 

In the sequel we put A(x,r) = {y £ R n ] \y — x\ < r} for any x g R n and 
r > 0. Let tu = inf{i > 0;Xf £ U c } for any open set U. And we suppose the 
operator © is elliptic on R n . 

Theorem 3. f(x) is a local bounded continuous function on R n . If for any 

x G R n , there exists a radius r(x), the mean value property holds: 

£ Q,t x [^ X tJ\= f( X ) here T * = T A(x,r(x))- (16) 

And r(x) is a measurable function of x and satisfies that for each x, there exists 
a bounded open set U X) x € U x and moreover r(y), y £ U x should satisfy the 
following two conditions: 

0<r{y)<dist(y,dU x ), (17) 

and 

inf{r(y);yeK}>0 (18) 

for all closed subsets K of U x with dist(K,dU x ) > 0. Then we can get: 

(i) . For each y £ U x the mean value property holds on the boundary: 

£l Ty [f{Xl)\ = f(y), here r y = inf{t > 0; X\ g U c x }. 

and furthormore: 

(ii) . f(x) is the viscosity solution of PDF (12). 

Proof, (i) =>■ (ii) is also based on the nonlinear Feynman-Kac formula for elliptic 
PDE(cf.[8]). So we sufficiently prove the first conclusion. 

For each y g U x , we define a sequence of stopping times for {Xf} by 
induction as follows: 

r =0 

r k - inf{* > r fe _ i; \Xf - Xy Tk i | > riXl J}; k > 1. 
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By the mean propertv (|16p . and the strong markovian property we can get 

a a x ^k i 

= £l Tk _ 1 [f(xy Tk _ 1 )}, 

then by induction we get 

^,J/TO1 = /(*)• 

In the following we will prove Tfc — > r y a.e. when k — > oo. Obviously 

so there exists a stopping time r s.t. Tk t T - If t 7^ t\,, then there exists e > 
s.t. 

dist(X^ , dU x ) > e for any k. 
Let rk = r(X*), according to the condition (fT5)) , we get there exists r > 0, 

r k > r fo r an y fe- 
lt means 

And since X t y is continuous, then Tk — > 00, which implies r y — 00. So 

P(r fc don't converge to r y ) < P(r a = 00). 
But for ([6]) is elliptic and U x is bounded, we have P(r y < 00) = 1. So 

P(t& converge to T y ) = 1. 

Then we get 

/(y) = ££ l7 J/TO] 

So we have finished the proof. □ 
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